Differential DNA methylation in iPSC-derived dopaminergic neurons: a step forward on the role of SNORD116 microdeletion in the pathophysiology of addictive behavior in Prader-Willi syndrome.

INTRODUCTION: A microdeletion including the SNORD116 gene (SNORD116 MD) has been shown to drive the Prader-Willi syndrome (PWS) features. PWS is a neurodevelopmental disorder clinically characterized by endocrine impairment, intellectual disability and psychiatric symptoms such as a lack of emotional regulation, impulsivity, and intense temper tantrums with outbursts. In addition, this syndrome is associated with a nutritional trajectory characterized by addiction-like behavior around food in adulthood. PWS is related to the genetic loss of expression of a minimal region that plays a potential role in epigenetic regulation. Nevertheless, the role of the SNORD116 MD in DNA methylation, as well as the impact of the oxytocin (OXT) on it, have never been investigated in human neurons. METHODS: We studied the methylation marks in induced pluripotent stem-derived dopaminergic neurons carrying a SNORD116 MD in comparison with those from an age-matched adult healthy control. We also performed identical neuron differentiation in the presence of OXT. We performed a genome-wide DNA methylation analysis from the iPSC-derived dopaminergic neurons by reduced-representation bisulfite sequencing. In addition, we performed RNA sequencing analysis in these iPSC-derived dopaminergic neurons differentiated with or without OXT. RESULTS: The analysis revealed that 153,826 cytosines were differentially methylated between SNORD116 MD neurons and control neurons. Among the differentially methylated genes, we determined a list of genes also differentially expressed. Enrichment analysis of this list encompassed the dopaminergic system with COMT and SLC6A3. COMT displayed hypermethylation and under-expression in SNORD116 MD, and SLC6A3 displayed hypomethylation and over-expression in SNORD116 MD. RT-qPCR confirmed significant over-expression of SLC6A3 in SNORD116 MD neurons. Moreover, the expression of this gene was significantly decreased in the case of OXT adjunction during the differentiation. CONCLUSION: SNORD116 MD dopaminergic neurons displayed differential methylation and expression in the COMT and SLC6A3 genes, which are related to dopaminergic clearance.

Single cell tracing of Pomc neurons reveals recruitment of 'Ghost' subtypes with atypical identity in a mouse model of obesity.

The hypothalamus contains a remarkable diversity of neurons that orchestrate behavioural and metabolic outputs in a highly plastic manner. Neuronal diversity is key to enabling hypothalamic functions and, according to the neuroscience dogma, it is predetermined during embryonic life. Here, by combining lineage tracing of hypothalamic pro-opiomelanocortin (Pomc) neurons with single-cell profiling approaches in adult male mice, we uncovered subpopulations of 'Ghost' neurons endowed with atypical molecular and functional identity. Compared to 'classical' Pomc neurons, Ghost neurons exhibit negligible Pomc expression and are 'invisible' to available neuroanatomical approaches and promoter-based reporter mice for studying Pomc biology. Ghost neuron numbers augment in diet-induced obese mice, independent of neurogenesis or cell death, but weight loss can reverse this shift. Our work challenges the notion of fixed, developmentally programmed neuronal identities in the mature hypothalamus and highlight the ability of specialised neurons to reversibly adapt their functional identity to adult-onset obesogenic stimuli.

Eomes-dependent mitochondrial regulation promotes survival of pathogenic CD4+ T cells during inflammation.

The mechanisms whereby Eomes controls tissue accumulation of T cells and strengthens inflammation remain ill-defined. Here, we show that Eomes deletion in antigen-specific CD4+ T cells is sufficient to protect against central nervous system (CNS) inflammation. While Eomes is dispensable for the initial priming of CD4+ T cells, it is required for long-term maintenance of CNS-infiltrating CD4+ T cells. We reveal that the impact of Eomes on effector CD4+ T cell longevity is associated with sustained expression of multiple genes involved in mitochondrial organization and functions. Accordingly, epigenetic studies demonstrate that Eomes supports mitochondrial function by direct binding to either metabolism-associated genes or mitochondrial transcriptional modulators. Besides, the significance of these findings was confirmed in CD4+ T cells from healthy donors and multiple sclerosis patients. Together, our data reveal a new mechanism by which Eomes promotes severity and chronicity of inflammation via the enhancement of CD4+ T cell mitochondrial functions and resistance to stress-induced cell death.

Macromolecular Complex Including MLL3, Carabin and Calcineurin Regulates Cardiac Remodeling.

BACKGROUND: Cardiac hypertrophy is an intermediate stage in the development of heart failure. The structural and functional processes occurring in cardiac hypertrophy include extensive gene reprogramming, which is dependent on epigenetic regulation and chromatin remodeling. However, the chromatin remodelers and their regulatory functions involved in the pathogenesis of cardiac hypertrophy are not well characterized. METHODS: Protein interaction was determined by immunoprecipitation assay in primary cardiomyocytes and mouse cardiac samples subjected or not to transverse aortic constriction for 1 week. Chromatin immunoprecipitation and DNA sequencing (ChIP-seq) experiments were performed on the chromatin of adult mouse cardiomyocytes. RESULTS: We report that the calcium-activated protein phosphatase CaN (calcineurin), its endogenous inhibitory protein carabin, the STK24 (STE20-like protein kinase 3), and the histone monomethyltransferase, MLL3 (mixed lineage leukemia 3) form altogether a macromolecular complex at the chromatin of cardiomyocytes. Under basal conditions, carabin prevents CaN activation while the serine/threonine kinase STK24 maintains MLL3 inactive via phosphorylation. After 1 week of transverse aortic constriction, both carabin and STK24 are released from the CaN-MLL3 complex leading to the activation of CaN, dephosphorylation of MLL3, and in turn, histone H3 lysine 4 monomethylation. Selective cardiac MLL3 knockdown mitigates hypertrophy, and chromatin immunoprecipitation and DNA sequencing analysis demonstrates that MLL3 is de novo recruited at the transcriptional start site of genes implicated in cardiomyopathy in stress conditions. We also show that CaN and MLL3 colocalize at chromatin and that CaN activates MLL3 histone methyl transferase activity at distal intergenic regions under hypertrophic conditions. CONCLUSIONS: Our study reveals an unsuspected epigenetic mechanism of CaN that directly regulates MLL3 histone methyl transferase activity to promote cardiac remodeling.

Single-cell RNA sequencing identifies senescence as therapeutic target in rhabdomyolysis-induced acute kidney injury.

BACKGROUND: The role of macrophages in the development of rhabdomyolysis-induced acute kidney injury (RM-AKI) has been established, but an in-depth understanding of the changes in the immune landscape could help to improve targeted strategies. Whereas senescence is usually associated with chronic kidney processes, we also wished to explore whether senescence could also occur in AKI and whether senolytics could act on immune cells. METHODS: Single-cell RNA sequencing was used in the murine glycerol-induced RM-AKI model to dissect the transcriptomic characteristics of CD45+ live cells sorted from kidneys 2 days after injury. Public datasets from murine AKI models were reanalysed to explore cellular senescence signature in tubular epithelial cells (TECs). A combination of senolytics (dasatinib and quercetin, DQ) was administered to mice exposed or not to RM-AKI. RESULTS: Unsupervised clustering of nearly 17 000 single-cell transcriptomes identified seven known immune cell clusters. Sub-clustering of the mononuclear phagocyte cells revealed nine distinct cell sub-populations differently modified with RM. One macrophage cluster was particularly interesting since it behaved as a critical node in a trajectory connecting one major histocompatibility complex class IIhigh (MHCIIhigh) cluster only present in Control to two MHCIIlow clusters only present in RM-AKI. This critical cluster expressed a senescence gene signature, that was very different from that of the TECs. Senolytic DQ treatment blocked the switch from a F4/80highCD11blow to F4/80lowCD11bhigh phenotype, which correlated with prolonged nephroprotection in RM-AKI. CONCLUSIONS: Single-cell RNA sequencing unmasked novel transitional macrophage subpopulation associated with RM-AKI characterized by the activation of cellular senescence processes. This work provides a proof-of-concept that senolytics nephroprotective effects may rely, at least in part, on subtle immune modulation.

The proteome and transcriptome of stress granules and P bodies during human T lymphocyte activation.

Stress granules (SGs) and processing bodies (PBs) are membraneless cytoplasmic assemblies regulating mRNAs under environmental stress such as viral infections, neurological disorders, or cancer. Upon antigen stimulation, T lymphocytes mediate their immune functions under regulatory mechanisms involving SGs and PBs. However, the impact of T cell activation on such complexes in terms of formation, constitution, and relationship remains unknown. Here, by combining proteomic, transcriptomic, and immunofluorescence approaches, we simultaneously characterized the SGs and PBs from primary human T lymphocytes pre and post stimulation. The identification of the proteomes and transcriptomes of SGs and PBs indicate an unanticipated molecular and functional complementarity. Notwithstanding, these granules keep distinct spatial organizations and abilities to interact with mRNAs. This comprehensive characterization of the RNP granule proteomic and transcriptomic landscapes provides a unique resource for future investigations on SGs and PBs in T lymphocytes.

Tamoxifen Accelerates Endothelial Healing by Targeting ERα in Smooth Muscle Cells.

RATIONALE: Tamoxifen prevents the recurrence of breast cancer and is also beneficial against bone demineralization and arterial diseases. It acts as an ER (estrogen receptor) α antagonist in ER-positive breast cancers, whereas it mimics the protective action of 17ß-estradiol in other tissues such as arteries. However, the mechanisms of these tissue-specific actions remain unclear. OBJECTIVE: Here, we tested whether tamoxifen is able to accelerate endothelial healing and analyzed the underlying mechanisms. METHODS AND RESULTS: Using 3 complementary mouse models of carotid artery injury, we demonstrated that both tamoxifen and estradiol accelerated endothelial healing, but only tamoxifen required the presence of the underlying medial smooth muscle cells. Chronic treatment with 17ß-estradiol and tamoxifen elicited differential gene expression profiles in the carotid artery. The use of transgenic mouse models targeting either whole ERα in a cell-specific manner or ERα subfunctions (membrane/extranuclear versus genomic/transcriptional) demonstrated that 17ß-estradiol-induced acceleration of endothelial healing is mediated by membrane ERα in endothelial cells, while the effect of tamoxifen is mediated by the nuclear actions of ERα in smooth muscle cells. CONCLUSIONS: Whereas tamoxifen acts as an antiestrogen and ERα antagonist in breast cancer but also on the membrane ERα of endothelial cells, it accelerates endothelial healing through activation of nuclear ERα in smooth muscle cells, inviting to revisit the mechanisms of action of selective modulation of ERα.

A novel locus on chromosome 1 underlies the evolution of a melanic plumage polymorphism in a wild songbird.

Understanding the mechanisms responsible for phenotypic diversification within and among species ultimately rests with linking naturally occurring mutations to functionally and ecologically significant traits. Colour polymorphisms are of great interest in this context because discrete colour patterns within a population are often controlled by just a few genes in a common environment. We investigated how and why phenotypic diversity arose and persists in the Zosterops borbonicus white-eye of Reunion (Mascarene archipelago), a colour polymorphic songbird in which all highland populations contain individuals belonging to either a brown or a grey plumage morph. Using extensive phenotypic and genomic data, we demonstrate that this melanin-based colour polymorphism is controlled by a single locus on chromosome 1 with two large-effect alleles, which was not previously described as affecting hair or feather colour. Differences between colour morphs appear to rely upon complex cis-regulatory variation that either prevents the synthesis of pheomelanin in grey feathers, or increases its production in brown ones. We used coalescent analyses to show that, from a 'brown' ancestral population, the dominant 'grey' allele spread quickly once it arose from a new mutation. Since colour morphs are always found in mixture, this implies that the selected allele does not go to fixation, but instead reaches an intermediate frequency, as would be expected under balancing selection.

PNPLA1 defects in patients with autosomal recessive congenital ichthyosis and KO mice sustain PNPLA1 irreplaceable function in epidermal omega-O-acylceramide synthesis and skin permeability barrier.

Autosomal recessive congenital ichthyosis (ARCI) is a heterogeneous group of monogenic genodermatoses that encompasses non-syndromic disorders of keratinization. The pathophysiology of ARCI has been linked to a disturbance in epidermal lipid metabolism that impaired the stratum corneum function, leading to permeability barrier defects. Functional characterization of some genes involved in ARCI contributed to the identification of molecular actors involved in epidermal lipid synthesis, transport or processing. Recently, PNPLA1 has been identified as a gene causing ARCI. While other members of PNPLA family are key elements in lipid metabolism, the function of PNPLA1 remained unclear. We identified 5 novel PNPLA1 mutations in ARCI patients, mainly localized in the putative active enzymatic domain of PNPLA1. To investigate Pnpla1 biological role, we analysed Pnpla1-deficient mice. KO mice died soon after birth from severe epidermal permeability defects. Pnpla1-deficient skin presented an important impairment in the composition and organization of the epidermal lipids. Quantification of epidermal ceramide species highlighted a blockade in the production of ω-O-acylceramides with a concomitant accumulation of their precursors in the KO. The virtually loss of ω-O-acylceramides in the stratum corneum was linked to a defective lipid coverage of the resistant pericellular shell encapsulating corneocytes, the so-called cornified envelope, and most probably disorganized the extracellular lipid matrix. Finally, these defects in ω-O-acylceramides synthesis and cornified envelope formation were also evidenced in the stratum corneum from PNPLA1-mutated patients. Overall, our data support that PNPLA1/Pnpla1 is a key player in the formation of ω-O-acylceramide, a crucial process for the epidermal permeability barrier function.

Lowering relative humidity level increases epidermal protein deimination and drives human filaggrin breakdown.

BACKGROUND: Deimination (also known as citrullination), the conversion of arginine in a protein to citrulline, is catalyzed by a family of enzymes called peptidylarginine deiminases (PADs). Three PADs are expressed in the epidermis, one of their targets being filaggrin. Filaggrin plays a central role in atopic dermatitis and is a key protein for the epidermal barrier. It aggregates keratins and is cross-linked to cornified envelopes. Following its deimination, it is totally degraded to release free amino acids, contributing to the natural moisturizing factor (NMF). The mechanisms controlling this multistep catabolism in human are unknown. OBJECTIVE: To test whether external humidity plays a role, and investigate the molecular mechanisms involved. METHODS: Specimens of reconstructed human epidermis (RHEs) produced in humid or dry conditions (>95% or 30-50% relative humidity) were compared. RESULTS: RHEs produced in the dry condition presented structural changes, including a thicker stratum corneum and a larger amount of keratohyalin granules. The transepidermal water loss and the stratum corneum pH were decreased whereas the quantity of NMF was greater. This highly suggested that filaggrin proteolysis was up-regulated. The expression/activity of the proteases involved in filaggrin breakdown did not increase while PAD1 expression and the deimination rate of proteins, including filaggrin, were drastically enhanced. Partial inhibition of PADs with Cl-amidine reversed the effect of dryness on filaggrin breakdown. CONCLUSION: These results demonstrate the importance of external humidity in the control of human filaggrin metabolism, and suggest that deimination plays a major role in this regulation.

SNP discovery and genetic mapping using genotyping by sequencing of whole genome genomic DNA from a pea RIL population.

BACKGROUND: Progress in genetics and breeding in pea still suffers from the limited availability of molecular resources. SNP markers that can be identified through affordable sequencing processes, without the need for prior genome reduction or a reference genome to assemble sequencing data would allow the discovery and genetic mapping of thousands of molecular markers. Such an approach could significantly speed up genetic studies and marker assisted breeding for non-model species. RESULTS: A total of 419,024 SNPs were discovered using HiSeq whole genome sequencing of four pea lines, followed by direct identification of SNP markers without assembly using the discoSnp tool. Subsequent filtering led to the identification of 131,850 highly designable SNPs, polymorphic between at least two of the four pea lines. A subset of 64,754 SNPs was called and genotyped by short read sequencing on a subpopulation of 48 RILs from the cross 'Baccara' x 'PI180693'. This data was used to construct a WGGBS-derived pea genetic map comprising 64,263 markers. This map is collinear with previous pea consensus maps and therefore with the Medicago truncatula genome. Sequencing of four additional pea lines showed that 33 % to 64 % of the mapped SNPs, depending on the pairs of lines considered, are polymorphic and can therefore be useful in other crosses. The subsequent genotyping of a subset of 1000 SNPs, chosen for their mapping positions using a KASP™ assay, showed that almost all generated SNPs are highly designable and that most (95 %) deliver highly qualitative genotyping results. Using rather low sequencing coverages in SNP discovery and in SNP inferring did not hinder the identification of hundreds of thousands of high quality SNPs. CONCLUSIONS: The development and optimization of appropriate tools in SNP discovery and genetic mapping have allowed us to make available a massive new genomic resource in pea. It will be useful for both fine mapping within chosen QTL confidence intervals and marker assisted breeding for important traits in pea improvement.

Shotgun assembly of the complete mitochondrial genome of the neotropical cracker butterfly Hamadryas epinome.

The complete mitochondrial genome of the cracker butterfly Hamadryas epinome (C. Felder and R. Felder, 1867) (Lepidoptera: Nymphalidae: Biblidinae) has been sequenced using a genome-skimming approach on an Illumina Hiseq 2000 platform. The mitochondrial genome of H. epinome was determined to be 15,207 bp long and presents an organization similar to other Ditrysia mitogenomes. A non-coding poly-AT region of uncertain length is present at position 6180.

Sequencing of the mitochondrial genome of the avocado lace bug Pseudacysta perseae (Heteroptera, Tingidae) using a genome skimming approach.

Lace bugs (Tingidae) are a family of phytophagous heteropterans, some of which are important agricultural and forestry pests. They currently comprise around 2500 species distributed worldwide, for which only one mitochondrial genome has been described so far. We sequenced the complete mitochondrial genome and the nuclear ribosomal gene segment of the avocado lace bug Pseudacysta perseae using a genome skimming approach on an Illumina Hiseq 2000 platform. Fifty-four additional heteropteran mitogenomes, including the one of the sycamore lace bug Corythucha ciliata, were retrieved to allow for comparisons and phylogenetic analyses. P. perseae mitochondrial genome was determined to be 15,850 bp long, and presented the typical organisation of insect mitogenomes. The phylogenetic analysis placed P. perseae as a sister to C. ciliata but did not confirm the monophyly of Miroidae including Tingidae. Our results contradicted widely accepted phylogenetic hypothesis, which highlights the limits of analyses based on mitochondrial data only. Shotgun sequencing approaches should provide substantial improvements in harmonizing mitochondrial and nuclear databases.

From museums to genomics: old herbarium specimens shed light on a C3 to C4 transition.

Collections of specimens held by natural history museums are invaluable material for biodiversity inventory and evolutionary studies, with specimens accumulated over 300 years readily available for sampling. Unfortunately, most museum specimens yield low-quality DNA. Recent advances in sequencing technologies, so called next-generation sequencing, are revolutionizing phylogenetic investigations at a deep level. Here, the Illumina technology (HiSeq) was used on herbarium specimens of Sartidia (subfamily Aristidoideae, Poaceae), a small African-Malagasy grass lineage (six species) characteristic of wooded savannas, which is the C3 sister group of Stipagrostis, an important C4 genus from Africa and SW Asia. Complete chloroplast and nuclear ribosomal sequences were assembled for two Sartidia species, one of which (S. perrieri) is only known from a single specimen collected in Madagascar 100 years ago. Partial sequences of a few single-copy genes encoding phosphoenolpyruvate carboxylases (ppc) and malic enzymes (nadpme) were also assembled. Based on these data, the phylogenetic position of Malagasy Sartidia in the subfamily Aristidoideae was investigated and the biogeographical history of this genus was analysed with full species sampling. The evolutionary history of two genes for C4 photosynthesis (ppc-aL1b and nadpme-IV) in the group was also investigated. The gene encoding the C4 phosphoenolpyruvate caroxylase of Stipagrostis is absent from S. dewinteri suggesting that it is not essential in C3 members of the group, which might have favoured its recruitment into a new metabolic pathway. Altogether, the inclusion of historical museum specimens in phylogenomic analyses of biodiversity opens new avenues for evolutionary studies.

Shotgun assembly of the assassin bug Brontostoma colossus mitochondrial genome (Heteroptera, Reduviidae).

The complete mitochondrial genome of the assassin bug Brontostoma colossus (Distant, 1902) (Heteroptera: Reduviidae) has been sequenced using a genome-skimming approach on an Illumina Hiseq 2000 platform. Fifty-four additional heteropteran mitogenomes, including five assassin bug species, were retrieved to allow for comparisons and phylogenetic analyses. The mitochondrial genome of B. colossus was determined to be 16,625 bp long, and consists of 13 protein-coding genes (PCGs), 23 transfer-RNA genes (tRNAs), two ribosomal-RNA genes (rRNAs), and one control region. The nucleotide composition is biased toward adenine and thymine (A+T=73.4%). Overall, architecture, nucleotide composition and genome asymmetry are similar among all available assassin bug mitogenomes. All PCGs have usual start-codons (Met and Ile). Three T and two TA incomplete termination codons were identified adjacent to tRNAs, which was consistent with the punctuation model for primary transcripts processing followed by 3' polyadenylation of mature mRNA. All tRNAs exhibit the classic clover-leaf secondary structure except for tRNASer(AGN) in which the DHU arm forms a simple loop. Two notable features are present in the B. colossus mitogenome: (i) a 131 bp duplicated unit including the complete tRNAArg gene, resulting in 23 potentially functional tRNAs in total, and (ii) a 857 bp duplicated region comprising 277 bp of the srRNA gene and 580 bp of the control region. A phylogenetic analysis based on 55 true bug mitogenomes confirmed that B. colossus belongs to Reduviidae, but contradicted a widely accepted hypothesis. This highlights the limits of phylogenetic analyses based on mitochondrial data only.

Gonad transcriptome analysis of pearl oyster Pinctada margaritifera: identification of potential sex differentiation and sex determining genes.

BACKGROUND: Black pearl farming is based on culture of the blacklip pearl oyster Pinctada margaritifera (Mollusca, lophotrochozoa), a protandrous hermaphrodite species. At first maturation, all individuals are males. The female sex appears progressively from two years old, which represents a limitation for broodstock conditioning for aquaculture production. In marine mollusks displaying hermaphroditic features, data on sexual determinism and differentiation, including the molecular sex determining cascade, are scarce. To increase genomic resources and identify the molecular mechanisms whereby gene expression may act in the sexual dimorphism of P. margaritifera, we performed gonad transcriptome analysis. RESULTS: The gonad transcriptome of P. margaritifera was sequenced from several gonadic samples of males and females at different development stages, using a Next-Generation-Sequencing method and RNAseq technology. After Illumina sequencing, assembly and annotation, we obtained 70,147 contigs of which 62.2% shared homologies with existing protein sequences, and 9% showed functional annotation with Gene Ontology terms. Differential expression analysis identified 1,993 differentially expressed contigs between the different categories of gonads. Clustering methods of samples revealed that the sex explained most of the variation in gonad gene expression. K-means clustering of differentially expressed contigs showed 815 and 574 contigs were more expressed in male and female gonads, respectively. The analysis of these contigs revealed the presence of known specific genes coding for proteins involved in sex determinism and/or differentiation, such as dmrt and fem-1 like for males, or foxl2 and vitellogenin for females. The specific gene expression profiles of pmarg-fem1-like, pmarg-dmrt and pmarg-foxl2 in different reproductive stages (undetermined, sexual inversion and regression) suggest that these three genes are potentially involved in the sperm-oocyte switch in P. margaritifera. CONCLUSIONS: The study provides a new transcriptomic tool to study reproduction in hermaphroditic marine mollusks. It identifies sex differentiation and potential sex determining genes in P. margaritifera, a protandrous hermaphrodite species.

Fast assembly of the mitochondrial genome of a plant parasitic nematode (Meloidogyne graminicola) using next generation sequencing.

Little is known about the variations of nematode mitogenomes (mtDNA). Sequencing a complete mtDNA using a PCR approach remains a challenge due to frequent genome reorganizations and low sequence similarities between divergent nematode lineages. Here, a genome skimming approach based on HiSeq sequencing (shotgun) was used to assemble de novo the first complete mtDNA sequence of a root-knot nematode (Meloidogyne graminicola). An AT-rich genome (84.3%) of 20,030 bp was obtained with a mean sequencing depth superior to 300. Thirty-six genes were identified with a semi-automated approach. A comparison with a gene map of the M. javanica mitochondrial genome indicates that the gene order is conserved within this nematode lineage. However, deep genome rearrangements were observed when comparing with other species of the superfamily Hoplolaimoidea. Repeat elements of 111 bp and 94 bp were found in a long non-coding region of 7.5 kb, as similarly reported in M. javanica and M. hapla. This study points out the power of next generation sequencing to produce complete mitochondrial genomes, even without a reference sequence, and possibly opening new avenues for species/race identification, phylogenetics and population genetics of nematodes.

Genome skimming by shotgun sequencing helps resolve the phylogeny of a pantropical tree family.

Whole genome sequencing is helping generate robust phylogenetic hypotheses for a range of taxonomic groups that were previously recalcitrant to classical molecular phylogenetic approaches. As a case study, we performed a shallow shotgun sequencing of eight species in the tropical tree family Chrysobalanaceae to retrieve large fragments of high-copy number DNA regions and test the potential of these regions for phylogeny reconstruction. We were able to assemble the nuclear ribosomal cluster (nrDNA), the complete plastid genome (ptDNA) and a large fraction of the mitochondrial genome (mtDNA) with approximately 1000×, 450× and 120× sequencing depth respectively. The phylogenetic tree obtained with ptDNA resolved five of the seven internal nodes. In contrast, the tree obtained with mtDNA and nrDNA data were largely unresolved. This study demonstrates that genome skimming is a cost-effective approach and shows potential in plant molecular systematics within Chrysobalanaceae and other under-studied groups.

Mass production of SNP markers in a nonmodel passerine bird through RAD sequencing and contig mapping to the zebra finch genome.

Here, we present an adaptation of restriction-site-associated DNA sequencing (RAD-seq) to the Illumina HiSeq2000 technology that we used to produce SNP markers in very large quantities at low cost per unit in the Réunion grey white-eye (Zosterops borbonicus), a nonmodel passerine bird species with no reference genome. We sequenced a set of six pools of 18-25 individuals using a single sequencing lane. This allowed us to build around 600 000 contigs, among which at least 386 000 could be mapped to the zebra finch (Taeniopygia guttata) genome. This yielded more than 80 000 SNPs that could be mapped unambiguously and are evenly distributed across the genome. Thus, our approach provides a good illustration of the high potential of paired-end RAD sequencing of pooled DNA samples combined with comparative assembly to the zebra finch genome to build large contigs and characterize vast numbers of informative SNPs in nonmodel passerine bird species in a very efficient and cost-effective way.

A DNA metabarcoding study of a primate dietary diversity and plasticity across its entire fragmented range.

In tropical regions, most primary ecosystems have been replaced by mosaic landscapes in which species must cope with a large shift in the distribution of their habitat and associated food resources. Primates are particularly vulnerable to habitat modifications. Most species persist in small fragments surrounded by complex human-mediated matrices whose structure and connectivity may strongly influence their dispersal and feeding behavior. Behavioral plasticity appears to be a crucial parameter governing the ability of organisms to exploit the resources offered by new matrix habitats and thus to persist in fragmented habitats. In this study, we were interested in the dietary plasticity of the golden-crowned sifaka (Propithecus tattersalli), an endangered species of lemur, found only in the Daraina region in north-eastern Madagascar. We used a DNA-based approach combining the barcoding concept and Illumina next-generation sequencing to (i) describe the species diet across its entire range and (ii) evaluate the influence of landscape heterogeneity on diet diversity and composition. Faeces from 96 individuals were sampled across the entire species range and their contents were analyzed using the trnL metabarcoding approach. In parallel, we built a large DNA reference database based on a checklist of the plant species of the Daraina region. Our results suggest that golden-crowned sifakas exhibit remarkable dietary diversity with at least 130 plant species belonging to 80 genera and 49 different families. We highlighted an influence of both habitat type and openness on diet composition suggesting a high flexibility of foraging strategies. Moreover, we observed the presence of numerous cultivated and naturalized plants in the faeces of groups living in forest edge areas. Overall, our findings support our initial expectation that P. tattersalli is able to cope with the current level of alteration of the landscape and confirm our previous results on the distribution and the dispersal ability of this species.